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We study the 56 Fe(> e , e") 56 Co cross section for the KARMEN neutrino spectrum. The Gamow- 
Teller contribution to the cross section is calculated within the shell model, while the forbidden 
transitions are evaluated within the continuum random phase approximation. We find a total cross 
section of 2.73 x 10 -40 cm 2 , in agreement with the data. 

PACS numbers: 26.50.+X, 23.40.-s, 21.60Cs, 21.60.Ka 
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The search for neutrino oscillations is one of the most 
promising experimental pathways for possible physics 
beyond the standard model of electroweak interac- 
tion. KARMEN is one of the pioneering low-energy 
accelerator-based experiments searching for these modes, 
ft consists of liquid scintillator housed in a 7000 ton 
shielding blockhouse (!]]. One of the experimental sig- 
nals which the KARMEN experiment can observe is the 
energy distribution of the electron emitted in charged- 
current reactions. This gives the experiment also the 
capability to obtain nuclear structure informations, and 
in fact the KARMEN collaboration has measured sev- 
eral (inclusive and exclusive) neutrino-induced charged 
and neutral current cross sections on 12 C 0, which is an 
important ingredient of the liquid scintillator. Recently 
the collaboration has also been able to determine inclu- 
sive (z/ e ,e~) cross sections for 13 C and 56 Fe. The latter 
might be quite important for two reasons: First, it pro- 
vides a benchmark to test the ability of nuclear models to 
describe neutrino-induced reactions on nuclei in the iron 
mass range as they are encountered in astrophysical sce- 
narios like the supernova collapse ||| . Second, future long 
baseline neutrino oscillation experiments like MINOS ||] 
or OMNIS f§ will use steel as detector material. Thus 
a rather reliable knowledge of the charged-current (and 
neutral-current) neutrino-induced cross section on 56 Fe 
is desirable to support the potential ability to observe 
supernova neutrinos by this detector. 

In this communication we present the first calculation 
(apart from a shell model estimate by Bugaev et al. |(|) of 
the 56 Fe(i/ e , e~) 56 Co cross section. The charged-current 
neutrino- induced cross section for a transition from a nu- 
clear state with angular momentum J and isospin quan- 
tum numbers T,Mt to a final nuclear state (J',T'Af^) 
is derived in Refs. Its evaluation reduces to the 

calculation of the various vector and axial vector multi- 
pole operators between appropriate nuclear many-body 
states. To derive these states we have adopted the follow- 
ing strategy, fn Ref. H it has been demonstrated that the 



inclusive cross section for forbidden transitions is rather 
well described within the continuum random phase ap- 
proximation (RPA). This model is, however, less reliable 
for the allowed transitions which in the long wavelength 
approximation are mediated by the Gamow- Teller oper- 
ator err. Here, the continuum RPA does not recover all 
nucleon-nucleon correlations necessary to reproduce the 
quenching of the Gamow- Teller strength. Such a task 
can be achieved in large-scale shell model calculations 
performed in the complete pf shell. On the other hand, 
at the higher neutrino energies present in the KARMEN 
i/ e -spectrum a replacement of the A = f + multipole op- 
erator by <jt is not quite appropriate. However, con- 
sidering the momentum-transfer dependence of the op- 
erator within shell model approaches represents a very 
formidable computational effort, even on modern paral- 
lel computers. To reliably estimate the contribution of 
the A = f + operator to the cross section, we have there- 
fore chosen a hybrid approach of shell model and RPA, 
as discussed below. 

Our shell-model calculation of the Gamow- Teller part 
to the 56 Fe(i^ e , e~) 56 Co cross section uses the Strassburg- 
Madrid codes ANTOfNE @ and NATHAN Q devel- 
oped by Caurier. The calculation has been performed 
within the complete pf shell adopting the recently de- 
veloped modified KB3 interaction. In this force slight 
monopole deficiencies encountered in the original KB3 
interaction have been corrected. Ref. |Q demonstrates 
that the modified KB3 version gives a fair account of 
nuclear spectra in the mass range A = 41 — 65 and 
reproduces the Gamow- Teller strength distributions in 
general very well, if the latter are scaled by the factor 
(0.74) 2 which accounts for the fact that complete Ohui 
shell-model calculations overestimate the experimental 
GT strength by a universal factor [f3]-[T^1 . We will adopt 
in the following the same universal scaling factor. With 
the code NATHAN it is now possible to calculate the to- 
tal GT strength for 56 Fe in the complete pf shell involving 
7 413488 J T = + configurations (this corresponds to an 
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m-scheme dimension of 501 million). We obtain a total 
GT_ strength (this is the direction in which a neutron is 
changed into a proton) of 8.85 units which agrees with the 
experimental value of 9.9 ± 2.4 |[6]]. However, a calcula- 
tion of the GT_ strength distribution Sgt(E), where E 
denotes the excitation energy in the final nucleus 56 Co, in 
the complete pf shell is still yet beyond present-day com- 
puter abilities. We have therefore truncated the model 
space in our calculation of the GT_ strength distribution 
of 56 Fe allowing a maximum of 5 particles to be excited 
from the / 2 orbital to the rest of the pf shell in the final 
nucleus. Furthermore the Ikeda sum rule has been ful- 
filled within the model space. The m-scheme dimension 
of this calculation is 1.6 million for 56 Fe and 19.8 million 
for 56 Co. In our truncated calculation we obtain a total 
GT_ strength of 9.3 units, showing that the calculation 
at this level of truncation is nearly converged. 

Due to the energy dependence of the phase space, the 
cross section is sensitive to the Gamow- Teller strength 
distribution rather than only to the total strength. 
Our shell-model calculation reproduces the experimen- 
tal GT_ strength distribution quite well, as can be seen 
in Fig. 2 of Rcf. Jl|. 

For the Gamow- Teller transition we then have the total 
cross section 

q2 cog 2 Q r 

a{E u ) = — ^ c J Pe E e F(Z + \,E e )S GT {E)dE e , 

(1) 

where E v is the neutrino energy, Gf is the Fermi con- 
stant, E e and p e are the energy and momentum of the 
outgoing electron, respectively, 8 C the Cabibbo angle, 
E = E u —E e due to energy conservation, and F(Z+1, E e ) 
the Fermi function which accounts for Coulomb distor- 
tion of the outgoing electron wave function. 

In the KARMEN experiment the neutrino beam has a 
Michel energy spectrum, 

96-E 2 

n(E„) = —f(M,-2E„) , (2) 

where M„ is the muon mass. Folding Eq. (1) with the 
spectrum (2) yields a partial 56 Fe(V e , e~) 56 Co cross sec- 
tion of 139xl0 -42 cm 2 for the Gamow- Teller transitions. 
As mentioned above, our treatment of the A = 1 + mul- 
tipole is only correct for small neutrino energies. To 
estimate the effect introduced by replacement of the 
A = 1 + operator by the momentum-transfer independent 
Gamow- Teller we have performed two calculations within 
the random phase approximation. In the first we adopted 
the complete A = 1 + multipole operator 0, in the second 
we replaced this operator by err. We then find partial 
56 Fe(^ e , e~) 56 Co cross sections of 341xl0~ 42 cm 2 and 
431 xlO -42 cm 2 , respectively. From these calculations 
we conclude that the replacement of the full A = 1 + op- 
erator by the Gamow- Teller operator increases the cross 



section by about 20%. If we correct our shell model esti- 
mate by this factor, we derive at our final result for the 
partial 56 Fe(^ e , e~) 56 Co cross section of HOxlO -42 cm 2 . 

The contributions of the other transitions to the cross 
section have been calculated using the continuum RPA. 
For nuclei with closed-shell configurations this model is 
well documented in the literature jlTj . We have recently 
generalized this formalism to allow also for partial occu- 
pancies in the parent ground state ||. From our shell- 
model calculation we find the following occupation num- 
bers for the pf shell orbitals: 117 / 2 = 5.30, n 3 / 2 = 0.32, 
n i/2 — 0.07 and 72.5/2 = 0.30 for protons n 7 / 2 = 7.29, 
n 3/2 — 1-57, nj/2 = 0.44 and n 5 / 2 = 0.70 for neu- 
trons. We assumed a completely occupied 40 Ca core. The 
single-particle energies were derived from an appropriate 
Woods-Saxon potential. We used the experimental val- 
ues for the proton and neutron separation energies. As 
the residual interaction we adopted the finite-range G- 
matrix derived from the Bonn potential |lq| . 

The other allowed transition is mediated by the + op- 
erator, which reduces to the Fermi transition in the limit 
of vanishing momentum transfer. The Fermi matrix ele- 
ment is easily obtained from isospin symmetry, yielding 
a total Fermi strength of Sj? = (N — Z) which is concen- 
trated in the isobaric analog state at 3.5 MeV in 56 Co. 
The Fermi contribution to the 56 Fe(^ e , e~) 56 Co cross sec- 
tion is then readily calculated within the RPA approach. 
We find a partial cross section of 53xl0 -42 cm 2 . 

The partial 56 Fe(i' e , e~) 56 Co cross sections for selected 
forbidden transitions are listed in Table |. From the for- 
bidden transitions only the J T = l - and 2~ dipole mul- 
tipoles are important. We find a summed cross section 
of 111 x 10 -42 cm 2 for the forbidden transitions, which 
is only slightly smaller than the contributions stemming 
from the allowed transitions. 

Adding up the allowed and forbidden transitions, we 
find a total 56 Fe(^ e , e") 56 Co cross section of 2.73 x 10" 40 
cm 2 , which agrees with the KARMEN result, 2.56 ± 
1.08(stat.) ± 0.43(syst.) x 10~ 40 cm 2 |§]. Fig. @ shows 
the differential 56 Fe(^ e , e~) 56 Co cross section as function 
of excitation energy in 56 Co. 

To test the importance of the partial occupancy for- 
malism for the forbidden transitions, we have repeated 
the continuum RPA calculation with a parent ground 
state whose occupation numbers have been taken from 
the independent-particle model. Thus, 717/2 = 6 for pro- 
tons, and 717/2 = 8 and 77.3/2 = 2 for neutrons. The 
respective partial cross sections are also listed in Table J. 
Confirming the results and arguments given in Ref. [BT, 
we observe that an improved description of the parent 
ground state, i.e. by considering partial occupancies, 
does not change the results noticeably. We now find a 
total cross section of 2.62 x 10 - 40 cm 2 . 

In summary, we have calculated the 56 Fe(V e , e~) 56 Co 
cross section for a neutrino spectrum corresponding to 
the KARMEN experiment. While the Fermi contribution 
for the transition to the isobaric analog state is trivial, 
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we calculate the Gamow- Teller piece to the cross section 
within a shell model approach which has been proven to 
correctly describe the GT_ strength distribution and to- 
tal strength for 56 Fe. For the forbidden transitions we 
employed the continuum RPA. Combining all relevant 
multipoles we find a total 56 Fe(^ e , e~) 56 Co cross section 
for the KARMEN neutrino spectrum of 2.73 x 1(T 40 cm 2 
which agrees with the experimental value obtained by the 
KARMEN collaboration, 256 ± 108(stat.) ± 43(syst.) x 
10~ 42 . Despite the rather large experimental uncertainty, 
this agreement is promising and shows that with the cur- 
rently chosen combination of nuclear models (shell model 
for allowed transitions and continuum RPA for forbid- 
den transitions) one should be able to reliably evaluate 
neutrino- induced cross sections on nuclei in the iron mass 
range as they are of interest in supernova simulations. 
Such calculations are in progress. 
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TABLE I. Partial cross sections o~j for the 56 Fe(u e , e~) 56 Co reaction induced by muon-decay-at-rest neutrinos. The second 
and third columns give the (continuum RPA) results calculated without and with the consideration of partial occupancies in 
the ground state. The Gamow- Teller result (J n = 1 + ) has been calculated within the shell model. The J = + multipole 
contribution reflects the Fermi transition to the isobaric analog state. The cross sections are in units of 10~ 42 cm 2 . 
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FIG. 1. Differential 56 Fe(> e 



3 Co cross section for the KARMEN neutrino spectrum as function of excitation energy in 



Co. The figure shows the allowed contributions, while the insert gives the contributions of the 1 and 2 multipolarities. 
The allowed contributions have been folded with a Gaussian of 0.5 MeV FHWM at energies below 5 MeV and with 1 MeV 
FHWM above 5 MeV. 
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